A Single Atom Transistor in a ID Optical Lattice 
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We propose a scheme utilising a quantum interference phenomenon to switch the transport of 
atoms in a ID optical lattice through a site containing an impurity atom. The impurity represents 
a qubit which in one spin state is transparent to the probe atoms, but in the other acts as a single 
atom mirror. This allows a single-shot quantum non-demolition measurement of the qubit spin. 

PACS numbers: 03.75.Lm, 42.50.-p, 03.67.Lx 



Coupling of a spin 1/2 system to Bosonic and 
Fermionic modes is one of the fundamental building 
blocks of quantum optics and solid state physics. Mo- 
tivated by the recent progress with cold atoms in ID 1], 
we consider a spin 1/2 atomic impurity which is used to 
switch the transport of either a ID Bose-Einstein Con- 
densate (BEC) or a ID degenerate Fermi gas initially 
situated to one side of the impurity. In one spin state 
the impurity is transparent to the probe atoms, whilst 
in the other it acts as single atom mirror, prohibiting 
transport via a quantum interference mechanism reminis- 
cent of electromagnetically induced transparency (EIT) 
2] (Fig. UJi). Observation of the atomic current pass- 
ing the impurity can then be used as a quantum non- 
demolition (QND) measurement [2j of its internal state, 
which can be seen to encode a qubit, \ip q ) = a\\) + /3||). 
If a macroscopic number of atoms pass the impurity, 
then the system will be in a macroscopic superposition, 

|*(t)) = a|T)|<Ar(*)) +/ 3 U)l^l( i ))> which can form the 
basis for a single shot readout of the qubit spin. Here, 
\<j) a {t)) denotes the state of the probe atoms after evolu- 
tion to time t, given that the qubit is in state a (Fig.[TJi). 
In view of the analogy between state amplification via 
this type of blocking mechanism and readout with single 
electron transistors (SET) used in solid state systems |4j , 
we refer to this setup as a Single Atom Transistor (SAT). 

We propose the implementation of a SAT using cold 
atoms in ID optical lattices [a, la Ll I3- We consider 
probe atoms b to be loaded in the lattice to the left of a 
site containing the impurity atom, which is trapped by 
a separate (e.g., spin-dependent JSJ) potential (Fig.^ 5 ). 
The passage of b atoms past the impurity q is then gov- 
erned by the spin-dependent effective collisional inter- 
action Hi n t = So- Uefi^blboqlqa. By making use of a 
quantum interference mechanism, we engineer complete 
complete blocking (effectively U c g — > oo) for one spin 
state and complete transmission (U e s — > 0) for the other. 
Below we first consider the detailed scattering processes 
involved in the transport of a single particle through 
the SAT, and then generalise this to interacting many- 
particle systems including a ID Tonks gas. 

The quantum interference mechanism needed to engi- 



(a) 



) \ / 



Wi 



FIG. 1: (a) A spin 1/2 impurity used as a switch: in one spin 
state it is transparent to the probe atoms, but in the other it 
acts as a single atom mirror, (b) Implementation of the SAT 
as a separately trapped impurity q with probe atoms b in an 
optical lattice. 



neer U e g can be produced using an optical or magnetic 
Feshbach resonance 9]. For the optical case a Raman 
laser drives a transition on the impurity site, 0, from the 
atomic state 6jgJ.|vac) via an off-resonant excited molec- 
ular state to a bound molecular state back in the lowest 



electronic manifold mj.|vac) (Fig. Et)- We denote the ef- 
fective two-photon Rabi frequency and detuning by Q a 
and A a respectively. The Hamiltonian for our system is 
then given (H = 1) by H — Hb + Hq, with 



H h = 



Ho = 



} j Qa (ml.q a bo + h.cj - A. a rh\rh a 
^\U qb , a blqlq a bo + U bm , a blmlm a b , (1) 



where the operators b obey the standard commutation 
(anti-commutation) relations for Bosons (Fermions). H^ 
gives a Hubbard Hamiltonian for the b atoms with 
tunnelling matrix elements J giving rise to a single 
Bloch band with dispersion relation e(k) = — 2J cos ka 
(a is the lattice spacing), and collisional interactions 
(which are non-zero only for Bosons) given by Ubb — 
4irh 2 ai,i„ J d 3 x |wj-(x)| /nib, where Wj(x) is the Wannier- 
function for a particle localized on site j, abb is the scat- 
tering length for b atoms and nib is their mass. Hq de- 
scribes the additional dynamics due to the impurity on 
site 0, where atoms b and q are converted to a molecular 
state with effective Rabi frequency Q a and detuning A a , 




FIG. 2: (a) The optical Feshbach setup couples the atomic 
state b ql\vac) (in a particular motional state quantised 
by the trap) to a molecular bound state of the Born- 
Oppenheimer potential, m£.|vac), with effective Rabi fre- 
quency Q a and detuning A CT . (b) A single atom passes 
the impurity (I — >III) via the two dressed states (II), +) = 
b ql-\vac} +mj.|vac) and | — ) = &o<jJ|vac) — m^vac). Quan- 
tum interference between the paths gives rise to an effective 
tunnelling rate J e s,a- 
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FIG. 3: (a) SAT transmission coefficients T = |/ (+) | 2 for a 
particle 6 as a function of its energy e(k) for Q/J = 4, A = 
0,U gb /J = (solid line), Q/J = 8, A/ J = 4,U qb /J = 2 
(dashed line), Q/J = 1,A = 0,U qb /J = 2 (dotted), and 
Q/J = 1,A = 0,U qb /J = (dash-dot), (b) The number 
of particles to the right of the impurity, Nn(t), from exact 
numerical calculations for Bosons in the limit U bb /J — » do 
(dashed lines) and Fermions (solid lines) in a ID Mott Insu- 
lator state with n — 1, for A = 0, Q/J — 0, 1, 2. 



and the last two terms describe background interactions, 
U a p,a for two particles a,/3 € {(ja-, &, wi}, which are typ- 
ically weak. This model is valid for U a p,J,Q,,A -C w, 
where u> is the energy separation between Bloch bands. 
Because the dynamics for the two spin channels q a can be 
treated independently, in the following we will consider 
a single spin channel, and drop the subscript a. 

For off-resonant laser driving (fi <§; |A|), the Fesh- 
bach resonance enhances the interaction between b and 
q atoms, giving the familiar result U e g = U q b + 2 /A. 
However, for resonant driving (A = 0) the physical mech- 
anism changes, and the effective tunnelling J c ff of an 
atom b past the impurity (Fig. |2b, I — > HI) is blocked 
by quantum interference. On the impurity site, laser 
driving mixes the states b q* \va,c) and 7rr |vac), forming 
two dressed states with energies e± = (U q b)/2± {U^ b /4 + 
C2 2 )V2 (Fig. 12b, II). The two resulting paths for a particle 
of energy e destructively interfere so that for large CI 3> J 
andC/ gb '=0, J cff = -J 2 /{e + tt)- J 2 /(e-ft) -> 0. This 
is analogous to the interference effect underlying EIT 
|2(, and is equivalent to having an effective interaction 
U e g — > co. In addition, if we choose A = —il 2 /U q b, the 
paths constructively interfere, screening the background 
interactions to produce perfect transmission (U e s — * 0). 

For a more detailed analysis, we solve the Lippmann- 
Schwinger equation exactly for scattering from the im- 
purity of an atom b with incident momentum k > in 
the lowest Bloch-band. The resulting forwards and back- 
wards scattering amplitudes, f^\k) respectively, are 



/ (±) (fc) 



v(k) 



±i 



(2) 



where the energy dependent interaction U c s = U q b + 
Q?/(e(k) — A) and the phase-velocity v(k) — de/dk = 
2Jasmka. The corresponding transmission probabili- 
ties, T(k) — \f( + \k)\ , are plotted in Fig. |3K as a func- 



tion oie(k) for various Q* and A. For intermediate il ~ J, 
these are Fano-profiles with complete reflection at e(k) = 
A and complete transmission at e(k) = A — il 2 /U q b- The 
SAT thus acts as an energy filter, which is widely tunable 
via the laser strength and detuning used in the optical 
Feshbach setup. For SI > 4J, T is approximately inde- 
pendent of k, and we recover the previous result, i.e., 
that transport can be completely blocked or permitted 
by appropriate selection of A. Note that this mecha- 
nism survives the generalisation to a multi-band Hubbard 
model, where higher energy Bloch bands are included in 
the calculation, and is resistant to loss process, which are 
discussed below. 

We now consider the full many-body dynamics of N 
probe atoms b initially prepared in the ground state in a 
trap (box) of M lattice sites on the left side of the impu- 
rity q. We are then interested in the expectation value 
of the steady state coherent current I = dNn/dt (where 
Nr = 2i>o bjbj is the number of particles on the right 
side of the impurity, see Fig- Et>), which depends on the 
laser parameters, the filling factor, n — N/M, and, for 
Bosons, the interaction strength, Ubb- We first consider 
the case of a dilute or noninteracting gas, before treating 
both interacting Bosons, and non-interacting Fermions 
with arbitrary n. 

For a dilute noninteracting Bose quasi-condensate 
(n <C I, Ubb = 0), or for any very dilute gas, (where 
the momentum distribution is very narrow), the be- 
haviour is very similar to that of a single particle. If 
the gas is quickly accelerated to a finite momentum k, 
e.g., by briefly tilting the lattice, then the atoms will co- 
herently tunnel through the impurity according to the 
scattering amplitudes f^(k). The resulting current 

I oc N |/^ + -'(fc)| v(k), where n is the initial filling fac- 
tor on the left of the impurity, and v{k) is the velocity of 
a Bloch- wave with momentum k. 
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FIG. 4: (a) The steady state current of b atoms through the 
impurity as function of the initial filling factor n for A = and 
Q/J = 0, 1, 2. The solid lines show the analytic result Jo for 
Fermions, whereas the dashed lines show the exact numerical 
result for hard-core Bosons with Ubb/J —* oo. For £2 = these 
results are indistinguishable, (b, c) The steady state current 
as function of the Rabi-frequency Q/J on resonance A = 
for (b) unit filling and (c) half filling. The solid lines show 
the analytic result for Fermions, whereas the dashed (dotted) 
lines give numerical results for Bosons with Ubb/J = 20 (4). 



For a Fermi gas the equations of motion are linear and 
may be solved exactly provided Ubm = U q b- Scattering 
from the impurity then occurs independently for each 
particle in the initial Fermi sea, with scattering ampli- 
tudes f^(k) for k < kp, where the Fermi momentum 
kp = nn. After a short transient period, on the order 
of the inverse tunnelling rate 1/J, the system establishes 
a roughly constant flux of particles through the impu- 
rity (Fig. |2t>) , with a time-averaged current for resonant 
driving A = given by 
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where Ep is the Fermi Energy in the initial state, /(e) is 
the density of states per site (on the left of the impurity) , 
2G\ = (l + £7 4 /4J 4 ) 1/2 ±l and V = e F /2J = sin 2 (mr/2). 

For a Tonks gas of strongly interacting Bosons 
(Ubb/J 3> 1 with n < 1) we expect to observe sim- 
ilar behaviour to that observed for Fermions. In this 
limit, double occupation of a site can be neglected, and 
the behaviour can be mapped onto Fermionic particles 
via a Jordan- Wigner transformation [lOj. The Hamil- 
tonian is then the same up to a nonlinear phase factor 
fi — > £2(— 1) B , which essentially causes £2 to change sign 
when a particle passes the impurity. The contribution 
of this phase factor should be small for weak coupling, 
£1 <C J, and also for strong coupling, where no particles 
will tunnel through the impurity, i.e., Nr ~ 0. 

For the general case of many Bosons we perform exact 
numerical integration of the time dependent Schrodinger 
equation for the Hamiltonian (Q using Vidal's algorithm 



for "slightly entangled quantum states" [llj . This algo- 
rithm selects adaptively a decimated Hilbert space on 
which a state is represented, by retaining at each time 
step only those basis states that carry the greatest weight 
in Schmidt decompositions taken from bipartite split- 
tings of the system between each two lattice sites. A 
sufficiently large decimated Hilbert space is then selected 
so that the results of the simulations are essentially exact. 
For each set of parameters we first prepared the initial 
state via an imaginary time evolution which found the 
ground state for atoms in a box trap on the left of the 
impurity. Then, considering initially a single impurity 
atom q on the site and unoccupied sites to the right 
of that site, we calculated the time evolution of the sys- 
tem until it had reached a quasi-steady state behaviour. 
These simulations allowed us to obtain the behavior for 
finite repulsion Ubb, and also to test the effects of the 
nonlinear phase factor (—1) Nr for strongly interacting 
Bosons, Ubb —> oo. 

In Fig. 0t> we plot the number of particles on the right 
of the impurity Nn(t) for Fermions and for Bosons with 
Ubb/J ~^ oo, starting from a Mott Insulator (MI) state 
with n = 1, for A = 0, Q/J = 0,1,2. For £1 = the 
results for Bosons and Fermions are identical, whilst for 
Q.J J = 1,2, we observe an initial period for the Bosons 
in which the current is similar to that for the Fermionic 
systems, after which the Bosons settle into a steady state 
with a significantly smaller current. The initial tran- 
sient period for the Bosons incorporates the settling to 
steady state of firstly the molecule dynamics, and sec- 
ondly the momentum distribution on the right of the im- 
purity. These transients are suppressed if O is ramped 
slowly to its final value from a large value Vt > 4 J. 

The dependence of the steady state current on the ini- 
tial filling factor n is depicted in Fig. QJi for resonant 
driving with Q/J = 0, 1, 2. For Q = U qb = U bm = 0, the 
current Iq = 2 J sin (mr/2)/Tr is identical for Fermions 
and hard-core Bosons (Ubb -» oo), as we expect from the 
exact correspondence given in this limit by the Jordan- 
Wigner transformation. For Fermions with weak, reso- 
nant laser driving, the main features of the Fano profile 
(Fig. |3t) are observed in correspondence with the inte- 
gral in (PI- For example, a plateau in Io(n) is observed 
near n = arccos(— A/2J)/tt = 1/2, as the Fermi Energy 
is raised past e ~ A = 0, which corresponds to the zero 
of the transmission probability T(e). Good agreement is 
also observed with the result for Bosons in this limit with 
n < 1/2, whilst for larger n Bosons are blocked better, 
with a factor of 2 — 3 in the steady state currents. 

The enhanced blocking for Bosons is also seen in 
Fig-Ek, where we plot the steady state current as a func- 
tion of O for resonant driving and n = 1. It is clear 
from these figures together that this difference is a fea- 
ture of the parameter regime n > 1/2, il ~ J, which 
can be linked directly to phase factor of (—l) Na that 
arises in the Jordan- Wigner transformation. As £7 is in- 



creased and fewer particles pass the impurity, the results 
for Fermions and Bosons again converge as expected. For 
small CI there arc small differences between Bosons with 
finite U^/ J — 4 and {U^/ J — > oo), with the currents al- 
ways lower than the equivalent fcrmionic current, owing 
largely to the smaller mean squared momentum in the 
initial state. For large driving, CI ^ 4J the basic inter- 
ference process is extremely efficient for both Bosons and 
Fermions, and we observe complete blocking or trans- 
mission of the current by quantum interference for the 
appropriate choice of A. 

In Fig. we investigate the time evolution of 30 hard- 
core Bosons (Ubb — > oo) in an initial MI state, which are 
released through a SAT which is switched at t = 0. For 
CI = 0, we see that as the gas expands the momentum dis- 
tribution becomes peaked as a quasi-condensate is formed 
with k = 7r/2a, which consists of a coherent superposi- 
tion of particles propagating to the right and holes prop- 
agating to the left as the MI state melts |f|. This mode 
grows outwards from the edge of the initial distribution, 
and contains at its peak ~ \/~N particles, as is expected 
for such dynamically formed quasi-condensates in a ID 
lattice 12] . In contrast, for CI/ J = 0.5, the momentum 
distribution is broader, and the quasi-condensate mode 
contains many fewer particles. The mode also consists 
of distinct branches, holes in the melting MI propagat- 
ing to the left and particles to the right, which are ini- 
tially coherent, but become decoupled at t ~ 12/ J. For 
larger CI this behaviour becomes more pronounced, and 
for CI > 4 J, the MI state essentially remains unchanged. 

The melting of a MI in this way can be used as the ba- 
sis for a convenient single-shot measurement of the spin 
state of q. If q is in a superposition of spin states, only one 
of which will permit transport of the b atoms, then after 
some propagation time, the system will be in a macro- 
scopic superposition of distinct quantum phases (MI and 
quasi-condensates). These are distinguishable because 
if the b atoms are released from the lattice, the quasi- 
condensate will produce an interference pattern, whereas 
the MI state will not. The visibility of the resulting pat- 
tern can thus be used to measure the qubit spin. 

A remarkable feature of the SAT is its resistance to 
both two- and three-body loss processes induced by the 
optical Feshbach resonance on the impurity site. Sponta- 
neous emissions from the off-resonant excited molecular 
state amount to a two-body loss process at a rate ~ 72s 
in the states |+) and | — ). These small rates are further 
suppressed in the blocking regime J, 72s <C CI, with the 
resulting decoherence rate 7d oc oc J 2 j2Bn/Cl 2 , with n the 
mean site occupation of the b atoms. Collisions of atoms 
b with molecules m [13j are strongly suppressed in the 
Tonks gas regime, as well as for Fermions. For a weakly 
interacting Bose gas the corresponding three-body loss 
rate, 73s, is again strongly suppressed in the blocking 
regime (J, 735 < CI) with 7 dcc oc J 4 7 3B n 2 /fi 4 . 

Parallels may be drawn between the SAT and other 
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FIG. 5: Exact numerical results showing the propagation of 
N = 30 hard-core Bosons in an initial MI state (n = 1) 
through a SAT with (a-c) CI = and (d-f) CI - 0.5J, with 
A = U q b — Ubm = 0. The plot shows (a,d) the momentum 
distribution, (b,e) the five largest eigenvalues, A m , of the sin- 
gle particle density matrix, (b\bj) and (c,f) the spatial density 
of the largest eigenmode (the quasi-condensate) as a function 
of time. Darker colours represent higher values. 



systems coupled to fermionic and bosonic modes. These 
include the QND-readout of a single photon in cavity- 
quantum electrodynamics [14|, electron counting statis- 
tics [13, and the transport of electrons past impurities 
such as quantum dots |l6| | (although there particles are 
normally initially present on both sides of the impurity). 
However, the long decoherence times for atoms in optical 
lattices imply coherent transport over longer timescales 
than is observed in these other systems, which are in- 
herently dissipative. In addition, blocking and/or energy 
filtering by one or more SATs could be applied as tools 
in the study of Bose and Fermi gases in a ID lattice. 
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